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ABSTRACT 

We searched for binary companions to 20 young brown dwarfs in the Upper 
Scorpius association (145 pc, 5 Myr, nearest OB association) with the the Laser 
Guide Star adaptive optics system and the facility infrared camera NIRC2 on 
the 10 m Keck II telescope. We discovered a 0.14" companion (20.9±0.4 AU) 
to the <0.1 M object SCH J16091837-20073523. From spectral deconvolution 
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of integrated-light near-IR spectroscopy of SCH1609 using the SpeX spectro- 
graph (Rayner et al. 2003), we estimate primary and secondary spectral types of 
M6±0.5 and M7±1.0, corresponding to masses of 79±17 Mj up and 55±25 Mj up 
at an age of 5 Myr and masses of 84±15 Mj up and 60±25 Mj up at an age of 
10 Myr. For our survey objects with spectral types later than M8, we find an 
upper limit on the binary fraction of <9% (1-cr) at separations of 10 - 500 AU. 
We combine the results of our survey with previous surveys of Upper Sco and 
similar young regions to set the strongest constraints to date on binary fraction 
for young substellar objects and very low mass stars. The binary fraction for low 
mass (<40 Mj up ) brown dwarfs in Upper Sco is similar to that for T dwarfs in the 
field; for higher mass brown dwarfs and very low mass stars, there is an excess of 
medium-separation (10-50 AU projected separation) young binaries with respect 
to the field. These medium separation binaries will likely survive to late ages. 

Subject headings: Upper Sco, brown dwarfs, planetary mass objects 



Introduction 



Numerous brown dwarf binaries have been discovered in the the field (Close et al.ll2003 



Burgasser et all 120031 : iBouv et all 120031 ; iBurgasser et al.l 120061 iLiu et all 120061 ) . Almost all 



of these have projected separations of <15 AU, with the majority having separations of <7 
AU. This tight binary distri bution was initial ly viewed as evidence for the ejection scenario 
of brown dwarf formation (IClose et al.l 120031 ) . In the ejection scenario, brown dwarfs are 
stellar embryos which are expelled from their natal subclusters due to interaction with other 
subcluster members, ther efore cutting off accretio n. Only tight brown dwarf binaries can 
survive an ejection event (IReipurth fc Clarkdl200ll ). 



In the last decade a population of wide (>15 AU separation) very low mass star, brown 
dwarf, and "plan etary mass" (<13 Mj„ p ) bir ~ 



nearby clusters (iLuhman 



Javawardhana fc Ivanovll2006l ; IClose et al. 
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see Table 1 for a list of all young <0.1 M binaries). These recent results 
suggest that the multiplicity properties of young (~few Myr) substellar objects in star- 
forming regions may be substantially different from the old (~few Gyr) field population. If 
common, these young binaries also provide serious constraints for current theories of brown 
dwarf f ormation, s ince such wide binaries cannot be formed by a non-dissipative ejection 
model ( jBatd 120091 ) . However, most of these objects were either discovered serendipitously, 
are from surveys with unpublished statistics, or are from surveys with very few objects of 
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comparable mass, so it is unknown how significant a population they form. Here, we conduct 
a systematic survey to search for such binaries in Upper Sco, the nearest OB association to 
the Earth. 



2. Sample Selection 



At an age of ~5 Myr and a distance of 145 pc (IPreibisch et al.ll2002[ ). the Upper Scorpius 
OB association is one of the nearest sites of ongoing star-formation and is intermediate in age 
between very young star- forming regions such as Taurus (<1 Myr) and somewhat older young 
field objects (~100 Myr). Additionally, Upper Sco is denser than nearby T associations such 
as Taurus and Chamaeleon but considerably less dense than high-mass star-forming regions 
such as the Trapezium in Orion (IPreibisch fc Mamaiekll2008l). Binarity of young objects may 



vary both as a function of age and environment (IPreibisch &: Mamajekl 120081 ). Certainly, 



the existence of yery young, wide binaries in <2 Myr star-formin g regions (e.g 



2004 



Todorov et al 



Lu 



nman 



AllerslboOjUjavawardhana fc Ivanovl bood : Idose et al J 120071 iKonopackvet alJl2007i 



2010h and the absence of such binaries in the field population (e.g 



Close et al. 



20031 ) suggests some evolution of brown dwarf binary properties must occur as a function 
of age. Thus, Upper Sco provides a key binarity data point, intermediate in both age and 
density. 

Some low m ass stars and high mass brown dwarfs in Upper Sco have already been stud- 
ied for binarity (IKraus et al.ll2005l 120081 ). Numerous binarity studies have been conducted 
which a re sensitive to very low substellar mass co mpanions for very young clusters such as 



(ILuhman et all 120051 ). NGC 1333 (IGreissl et al. 



aurus (jKraus et al.l 20061 ; iKonopackv et al.ll2007), Ch amaeleon (lAhmic et al.ll2007), IC 348 



20071 ). as well as Upper Sco (IKraus et al. 



2OO5J) - however, the least massive primaries observed in these surveys have generally been 
limited to higher mass brown dwarfs (>40 Mj up ). To date, only 18 objects with estimated 
masses <40 Mj up possess AO or space based observations for binarity which are published 
in surveys with well-defined contrast limits (7 objects from Kraus et al. 2005, 11 objects 
from Luhman et al. 2005). 

Here, we extend binarity results to lower mass brown dwarfs and planetary mass objects 
in Upper Sco. We surveyed a sample of 20 substellar objects in Upper Sco with reported 



spectral type s of M7.5 or 



confirmation ( Lodieu et al 



ater. These objects were selected from those with spectroscopic 



200 8J) fr o m the near-IR photometric and pro per motion su r veys o f 



Lodieu et all (120071 ): ISlesnick et all (120081 ). According to the models of iBaraffe et al.l (120031 ). 

these objects are all substellar. Indeed, these are the least massive objects currently known 
in Upper Sco, with estimated masses of <40 Mj up , thus this survey doubles the number of 
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young, low mass brown dwarfs imaged to search for binarity. At the 5 Myr age of Upper 
Sco, these objects are quite hot, hence their late M and early L spectral types. Eventually, 
these objects will c ool to become T dwarfs. We focus in particular on 18 objects selected 
from lLodieu et al.l (120081 ) which form a consistently selected and analyzed sample. Survey 
objects are listed in Table [2j 



Observations and Data Reduction 



We observed 20 objects with the facility i nfrared camera NIRC2 and the Laser Guide 
Star adaptive optics system fjBouy et al.ll2004t IWizinowich et al.ll2004l ) on the 10 m Keck II 
telescope. Observations were conducted on the nights of 2007-07-17, 2008-07-27, 2009-05-29, 
2009-05-30, and 2009-06-30. Conditions varied considerably between nights. We used the 
NIRC2 narrow camera, with a 9.963±0.005 mas pixel" 1 platescale and a 10.2"xl0.2" field 
of view. Search observations were conducted in the Ks filter (X ce ntrai = 2.146 /im). Objects 
were observed using a 3 point dither pattern, with a dither of 1-2.5" between positions. 
Observations are detailed in Table |3j The data were reduced in real time at the telescope 
using a custom IRAF pipeline. In cases where a candidate companion was detected at the 
telescope, immediate followup observations in J and H bands were then conducted. The 
observed object FWHM (Kg band) varied from 55 - 130 mas, with Strehl ratios in Ks 
varying from 6 to 29%. FWHMs and Strehl ratios were calculated using the standard Keck 
LGS routine nirc2strehl.pro. Most objects appeared slightly elongated in the direction of the 
tip-tilt reference star. We used a custom IDL pipeline for a final reduction of the data. The 
IDL pipeline corrects for on-chip distortion, flat fields, sky subtracts, and registers images 
using a cross-correlation algorithm. 



4. Candidate Selection Technique and Tentative Companion Candidates 

Images were visually inspected for candidate companions. A number of faint candidate 
companions to several survey objects were identified at separations of >1". To be consid- 
ered true companions, candidates must possess red colors similar to their primary and have 
common proper motion. Candidates to USco J160603. 75-221930.0 and USco J160723.82- 
221102.0 were reobserved 1 year after the initial observations and found to be background 
(i.e. not common proper motion objects). Colors for other candidates were checked in the 
ZYJHK bands of the UK Infrared Deep Sky Survey (UKIDSS) or in the Digital Sky Survey 
(DSS). The UKIDSS project is defined in Lawrence et al (2007). UKIDSS uses the UKIRT 
Wide Field Camera (WFCAM; Casali et al, 2007). The photometric system is described in 
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Hewett et al (2006), and the calibration is described in Hodgkin et al. (2009). The pipeline 
processing and science archive are described in Irwin et al (2009, in prep) and Hambly et al 
(2008). One candidate companion to USco J160830. 49-233511.0 was too faint to be detected 
in the UKIDSS data (K s ~19). This candidate will be reobserved at Keck in Spring 2011; 
however, given its faintness and wide separation (~5.4", ~780 AU), it is most likely back- 
ground. All of the other >l"candidates were detected with S/N > 10 in our Keck LGS AO 
data and were well-detected in the UKIDSS data. All UKIDSS detected objects were found 
to have colors significantly bluer than their primary; this is a clear sign that these objects 
are blue background objects as opposed to a red brown dwarf or planetary mass companion. 

Since our survey objects are quite faint and our AO correction is in general moderate, 
we did not attempt PSF subtraction to search for faint companions within 1" of the object. 
Most of our objects show some elongation towards the tip-tilt star. Additionally, image 
quality varied considerably between nights and during individual nights. Thus, it was not 
possible to build a reliable synthetic PSF for PSF subtraction. We note that our brighter 
targets had a number of superspeckles evident within 0.5" of the primary which can mimic 
the appearance of a companion. However, these superspeckles modulate with wavelength 
and also evolve as a function of time. By comparing multiple images taken at different times 
or wavelengths, it is almost always possible to distinguish speckles from real companions. In 
fact, we did initially flag a number of close-in candidates which proved to be speckles upon 
further examination. 



5. Discovery of a Brown Dwarf Companion to SCH 16091837-20073523 

A close candidate companion (0.14") was detected around SCH J16091837-20073523 
(henceforth SCH 1609-2009) with colors consistent with a young substellar object. Photom- 
etry and astrometry for this object is presented in Table HI Photometry and astrometry were 
determined using two different methods: 1) DAOPHOT psf-fitting photometry using IRAF 
and 2) synthetic psf-building photometry using BINFIT and StarFinder in IDL. 

For the DAOPHOT psf-fitting photometry, a background object in the field (with sepa- 
ration >5" from the primary and blue colors as expected for a background object) was used 
as a PSF for the allstar task. The PSF object was well detected in both H and K bands. 
However, our reduced AO correction in the blue J band relative to the H and K resulted in 
a lower signal to noise detection of the PSF star in the J band, diminishing our photometric 
accuracy in J band. 

Since we were somewhat concerned that the background object used as a PSF at >5" 
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might be affected by anisoplanetism, we also determined the the binary separ ation, position 



angle (PA), and flux ratio using the StarFinder PSF-fitting software package (jDiolaiti et al. 



20001 ) as an independent confirmation of the DAOPHOT results. StarFinder simultaneously 
solves for an empirical PSF model and the positions and fluxes of the binary components. Our 
J-band images had significantly poorer FWHMs and Strehl ratios such that StarFinder did 
not converge on a soluti on. Thus, for J ban d we instead used a three-component Gaussian 
model, as described by iDupuy et al.l (120091 ). to derive the binary parameters from PSF- 
fitting. The uncertainties were determined from the rms of the best-fit parameters for ou r 
individual dithered images. We adopted the astrometric calibration of iGhez et al.l ( 120081 ). 
with a pixel scale of 9.963±0.005 mas pixel -1 and an orientation for the detector's +|/-axis 
of +0?13±0?02 east of north. We applied the distortion correction developed by B. Cameron 
(2007, private communication), which changed our astrometry below the la level. Both 
DAOPHOT and StarFinder methods yielded nearly identical values (within the cited errors) 
for photometry and astrometry. 

We show the primary and companion on color-magnitude and color-color diagrams in 
Fig [21 The companion possesses very similar colors to its primary, suggesting that it is a 
true substellar companion. We estimated the likelihood that this companion is a background 
object using source counts from the 2MASS survey. Within 1 degree of the primary, 2MASS 
detects 527 objects with J of 13 mag or brighter, 506 objects with H of 12.4 or b r ighter , 
and 427 with K of 12 or brighter. Thus, adopting the approach of iBrandner et al.l (120001 ). 
in particular, their equation 1, we estimate the probability of finding an unrelated source at 
least as bright as the observed companion within 0.14" of the primary to be ~2.6xl0 -6 . 

SCH 1609-2007 was reobserved with NIRC2 at Keck II on 1 May 2010. The overall 
quality of the dataset was poor; however, we acquired sufficient data to demonstrate that the 
companion likely has common proper motion with the primary. Measuring centroid positions 
of the primary and companion (as the 2nd epoch data were not high enough quality for psf- 
fitting photometry), the companion moved by <0.7 pixels relative to the primary between 
epochs, consistent with the errors in our simple center-of-light centroiding. As no directly 
measured proper motion is av ailable for SCH 1609-2007 we adopt the mean value of (-11, -25) 
mas yr _1 for Upper Sco here (jde Bruijne et al.lll997l ; iPreibisch et al.lll998l ). At a distance of 
145 pc, parallax motion for Upper Sco is quite small - ~7 mas. As the parallax factor in 
the 2nd epoch observation was similar to that in the first, we neglect parallax here. With 
the ~10 mas pixel scale of narrow camera, we would have expected the companion to move 
~2.3 pixels relative to the primary between epochs if it was a background object at a much 
larger distance. Thus, this is likely a proper motion pair. 
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5.1. Spectroscopy and Spectral Type Estimates 



It has been noted that objects from ISlesnick et al. (120081) are considerably brighter 
than objects of similar spectral type from Lodieu et al. ( 20081 ). In fact, in some cases, the 
discrepancy is as much as 2 or 3 magnitudes, e.g. the M8 objects SCH J1622-1951 and 
USco J155419. 99-213543.1 in the sample for this survey . One possible explanation for this 
discrepancy is that the later type ISlesnick et al.l (120081 ) objects consist primarily of nearly 
equal mass binaries, such as SCH J1609 and likely SCH 1622-1951 as well. However, even 
after accounting for bin arity, SCH J1609-2 007 is still 2-3 mag overluminous compared to 
similar objects from the lLodieu et al.l ( 120081 ) sample. 



The discrepancy may also be due to systematic differences between optical and infrared 
spe ctral types for these o bjects, which are right at the M to L type spect ral transition . All o f 
the ISlesnick et al.l (120081 ) sources have optical spectral types while the lLodieu et al.l (120081 ) 
sources have infrared ones, so in effect we may be comparing apples vs. oranges. Thus, 
to further constrain the near-IR spectral type (and hence mass) of SCH J1609-2007AB we 
obtained integrated-light near- IR spectroscopy of SCH 1609-2007 on 2010 September 14 (UT) 
using the SpeX spectrograph ( Rayner et al. 20031 ) on the NASA Infrared Telescope Facility. 
A series of 12 exposures of 30 seconds each were taken, nodding along the slit, for a total 
integration time of six minutes. Our observations were taken at an airmass of 1.57, and the 
seeing recorded by the IRTF was / /9 The data were taken using the Low-Res prism with the 
CK/5 slit aligned with the parallactic angle, producing a 0.8-2.5 /im spectrum with a resolution 
(R=A/AA) of ~150. For telluric correction of our SCH1609-2007 spectrum, we observed a 
nearby A0V star, HD 149827 and obtained calibration frames (flats and arcs). The spectra 



were reduced using the facility reduction pipeline, Spextool (jCushing et al. 



200 4]), which 



i nclud es a correction for telluric absorption following the method described in IVacca et al. 
(120031 ). Spectra and spectral fits are presented in Fig [31 



SCH1609-2007AB was assigned a composite optical spectral type of M7.5 by lSlesnick et al 



(120081 ). We determine spectral types for each component by comparing our integrated light 
near-IR spectrum of SCH1609-2007AB to synthetic composites generated from template 
near-IR spectra of known members of Upper Scorpius (also taken with SpeX at the IRTF, at 
the same resolution as the SCH1609-2007AB spectrum, Brendan Bowler, priv comm). Our 
Upper Scorpius templates have optical spectral types ranging from M4 to M8.5. We verified 
that our templates have near-IR spectral types (calculated using the H 2 index of Allers et 
al. 2007) that agree to within 1 subtype with their optical types. 

To create our synthetic composite spectra, we first calculated synthetic photometry for 
each template using the J, H, and Ks filter profiles for NIRC2 and normalized each spectrum 
by the photometric flux density. We interpolated the templates to the same wavelength 
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grid as our spectrum of SCH1609-2007, and summed pairs of template spectra together, 
multiplying the later spectral type template by the flux ratio of the binary. Following the 
technique described in Cushing et al. (2008), we determined a multiplicative constant for 
each template in each band (J, H, and Ks), and computed a reduced y 2 performed over the 
wavelength ranges A = 0.96-1.3 fim, 1.48-1.8 /xm, and 2.05— 2.4/im. The best fitting template 
is the composite spectrum of UScoCTIO 75 (M6, Ardila et al. 2000; Preibisch et al. 2002) 
and DENIS-P J155605.0-210646 (M7; Martin et al. 2004, Slesnick et al. 2008). We assign 
spectral types of M6.0 ± 0.5 to SCH1609-2007A and M7.0 ±1.0 to SCH1609-2007B, where 
uncertainties are determined from the spectral types of synthetic composite spectra where 



Y 2 > X min + 1' significantly earlier than the combined M7.5 spectral type from lSlesnick et al. 
(r2008l ). 



5.2. Mass Estimates 



We estimate t he masses and effectiv e temperatures of SCH 1609- 2007 AB based on the 
DUSTY models of Chabrier et al. (hood ) and the temperature scale of lLuhmanl (120041) . The 
age of Upper Scorpius has been measured as 5 Myr, with a spread of up to 2 (IPreibisch Sz Zinnecker 
19991 ; ISlesnick et al.l 120081 ). but more recent work suggests ages as old as 10 Myr (Eric Ma- 
majek, private communication). 

Thus to account for age spread, we est imate masses at discrete ages of 5 and 10 Myr using 
the DUSTY models (jChabrier et al.l 120001) and a t an a ge range of 5±1 Myr using dust-free 
models from the same group faaraffe et al1ll998l . \200M For single age mass estimates, we 
simulated the spectral type range of each object with an input distribution of 10 6 Gaussian- 
distributed spectral type values centered on the measured spectral type and with a set to the 
error on the measured spectral type. We the n convert e d this input distribution to effective 
te mperatures using the t emperature scale of iLuhmanl ( 120041 ) and to estimated mass using 
the I Chabrier et al.l (120001 ) models. The estimated mass of each object was set to the mean 
of the output distribution and the error on the mass was set to the standard deviation of the 
output distribution. Via this method, we estimate primary and secondary masses of 79±17 



1 While these two sets of models differ in colors due to different atmospheric compositions (dust grains or 
the lack thereo f), they produce th e same bolometric luminosities and effective temperatures as a function of 
age and mass (jBaraffe et al.l 120021) . In particular, since isochrones are only defined at 1, 5 , and 10 Myr for 
the DUSTY models as opposed to a much denser grid of isochrones for the iBaraffe et al.l (|l998j ) models - 
and the authors sug gest caution using isochrones with ages <1 Myr, we have chosen to interpolate from the 



Baraffe et al.l (|1998l ) models when deriving masses at a range of ages to avoid inaccuracies from interpolating 



from I Myr isochrones. 
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Mj up and 55±25 Mj up at an age of 5 Myr and masses of 84±15 Mj up and 60±25 Mj up at an 
age of 10 Myr. 

For mass estimates for a range of ages, we simulated input spectral type and ages with 
an input distribution of 3xl0 4 Gaussian-distributed spectral type and age values. As before, 
the spectral type values were centered on the measured spectral type, with a set to the 
error on the measured spectral type. The center of th e age distribution was set as 5 Myr, 



with q=l M yr. Interpolation with the iLuhmanl (120041 ) temperature scale and iBaraffe et al. 



(119981 . 120021 ) models was performed to convert from spectral type to estimated mass for the 
distribution. Since the output distribution is somewhat asymmetric, we adopt the median 
as the best mass estimate and again use the standard deviation to set the error. Via this 
method, we estimate primary and secondary masses of 79±21 Mj up and 60±31 Mj up . Thus, 
for the range of ages that are realistic for this binary, the uncertainty in the measured spectral 
type dominates the mass estimate above and beyond any uncertainty in the age. 



5.3. Orbital Period Estimates 

We estimate the semimajor axis of SCH 1609-2007AB's orbit from its observed separa- 
tion. Assuming a uniform eccentricity distribution between < e < 1 and random viewing 
angles, Dupuy & Liu (2010) compute a median correction factor between projected sepa- 
ration and semimajor axis of l.lOig^g (68.3% confidence limits). Using this, we derive a 
semimajor axis of 23.0~^\ 9 5 ° AU for SCH 1609AB based on the observed separation in June 
2009. These correspond to an orbital period estimate of 310^2ii y r ; for an assumed total 
system mass of 134±30 Mj up . 



6. Achieved Contrasts and Limits on Minimum Detectable Companion Masses 



The 5-cr contrast curves for our core sample of 18 objects from lLodieu et al.l (120081 ) are 
presented in Figure |U Noise levels after data reduction were calculated as a function of 
radius by calculating the standard deviation in an annulus (with width equal to the FWHM 
of the PSF) centered on that radius. Noise curves were then converted to contrast in Amag 
by dividing by the measured peak pixel value of the object. Cont r asts w ere converted 
into absolute magnitudes using photometry reported in lLodieu et al.l (120081 ) and adopting 
a distance of 145 pc f or Upper Sco. A filt er transform was calculated from K to K s band 
using the spectra from lLodieu et al.l (120081 ) . Absolute magnitudes of the faintest detectable 
objects are also presented in Figure HI A table of contrast values at separations of 0.07, 0.2, 
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and 0.5"is presented in Table |5j 

To test the fidelity of our contrast curves, we inserted and retrieved simulated objects 
in our data. Objects were simulated as 2 dimensional gaussians with FWHMs from fits to 
the primary using the IDL routine GAUSSFIT2D and contrasts from our measured contrast 
curves. Objects simulated with contrasts from our measured curves were retrieved with 
S/N>5 for all survey targets down to separations of 0.07". For separations down to 0.06", 
simulated objects were retrieved for half of our survey targets. No simulated objects were 
retrieved at separations < 0.05". Thus, we conclude that our measured contrast curves are 
a reliable estimate of the detectable contrasts for potential companions down to separations 
of 0.07". 

We note that these contrasts do not take into account confusion between potential com- 
panions and speckles. Our brighter targets had a number of superspeckles evident within 
0.5" of the primary which can mimic the appearance of a companion. However, these su- 
perspeckles modulate with wavelength and also evolve as a function of time. By comparing 
multiple images taken at different times or wavelengths, it is almost always possible to dis- 
tinguish speckles from real companions. Thus, since we can distinguish between the two, 
we believe that our contrast curves adequately measure obtained contrasts for this survey, 
despite potential speckle confusion. 



Contrasts were converted to minimum detectable mass ratios using the models of lChabrier et al 



(120001 ) at an adopted age of 5 Myr and assuming a similar bolometric correction (i.e. a similar 
spectral types for both objects) between each target and any potential companion (Figure [5] 
and Table [5j). We note that for the best 75% of our data we are complete for all binaries 
with q>0.8 at separations >10 AU and all binaries with q>0.2 at separations >50 AU. 



7. Discussion 

7.1. Measured Binary Fraction 

We note that our newly discovered binary, SCH 1609 AB, is consistent with other young, 
wide very low mass binaries discovered, with a wide (>10 AU) separation and nearly equal 
mass ratio (q~0.7). With only one companion detected as part of our survey, we cannot 
place any new constraints on the mass ratio distribution or separation distribution for young 
brown dwarf companions. However, we have surveyed the largest sample to date of young 
brown dwarfs with estimated masses <40 Mj up and can strengthen constraints on the binary 
fraction (10 - 500 AU) of young objects in this mass range. We find an upper limit on the 
binary fraction (10 - 500 AU) of 9% (1-<t) for the 18 objects we surveyed from Lodieu et al. 
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2008 (calculated via the method of Burgasser et al. 2003). (We exclude the sources observed 
from the Slesnick et al. 2008 sample since they appear so much brighter than the Lodieu et 
al. sources and likely have masses > 40 Mj up ). 



7.2. Methods for Statistical Comparisons between Samples 

Given a sample of objects with true binary fraction eun-, the probability density of finding 
k binaries among n objects observed is given by the binomial distribution: 



f(k; n, e bm ) = — — — e£ n (l - e bin ) n ~ k (1) 

In our case, we would like to invert this probability density in order to obtain the 
probability of the sample having a given binary fraction e^ n in the case that we measure k 
binaries among n objects. To estimate the true binary fraction for our sample, we can then 
derive a confidence interval (presented as I— a intervals here) around the maximum of this 
probability density in which we expect the true binary fraction to reside. The probability 
density function and the resulting confidence intervals can either be calculated numerically 
(via e.g. the method of Burgasser et al. 2003) or by Bayesian posterior inference (see e.g. 
Sivia et al. 2006 and Cameron 2010). 

Quantitatively comparing the binary fractions (with error bars included from confidence 
intervals) from sample to sample requires some additional mathematical machinery. In some 
cases the confidence intervals overlap for binary fractions derived for different samples - 
however, it is not immediately clear how statistically significant this correlation is. In com- 
paring two samples of objects the question we wish to answer is: are they drawn from the 
same binomial distribution with e^n or from different distributions? To determine the likeli- 
hood that two samples are drawn fro m the same b i nomia l distribution, we adopted both the 



Fischer exact te s t met hod (used by lAhmic et al.l (120071 ) and described in the appendix of 



Brandeker et al.l ( 120061 )) as well as a Bayesian approach, derived below (derivation adopted 



from Carpenter 2009). 

According to Bayes' Theorem: 

probihypothesis; data, I) oc prob(data; hypothesis, I) x probihypothesis; I) (2) 

where (in this case) "I" is prior information, "data" is our measured sample, and "hy- 
pothesis" is the hypothesis (e.g. in this case we hypothesize that for brown dwarfs, phe- 
nomenon of binarity can be modeled as a binomial distribution with binary probability e&j n ). 
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In Bayesian terms, prob(hypothesis;I) is the prior probability and represents what we ini- 
tially know regarding the truth or falseness of the hypothesis while prob(data;hypothesis,I) 
is the likelihood function and gives the likelihood of each possible experimental outcome 
given the adopted model for the data. Combining the two gives prob(hypothesis;data,I), the 
posterior probability - the likelihood of a given model, in light of the measured data. 

In this case, we would like to derive the posterior probability density not for each 
individual sample but for the difference of the two: 



8 — tbinl — £bin2 (3) 

To do this, we first must derive the posterior probability distributions appropriate for 
each of the two binomial distributions we are comparing. Our likelihood function is again 
given by the binomial distribution, where eu n is the true binary fraction, for the sample, n 
is the total number of objects observed, and k is the number of binaries found: 



prob( data; hypothesis, I) = f(k;n,p) = — rrrP fc (l ~ P) n k (4) 

k'.yn — k)\ 

For the prior probability, prob(hypothesis;I), we simply adopt a uniform distribution 
from to 1, i.e. the binary fraction must be between and 1. This can also be written 
in terms of the beta distribution, a special case of the Dirichlet distribution with only two 
parameters defined on the interval (0,1): 

Adopting a = (3 =1, f(x;l,l) reduces to a uniform distribution, thus: 



probihypothesis; I) = 1 in the interval 0, 1 = Beta(l, 1) (6) 

The advantage of choosing the Beta distribution to represent the prior probability is that 
the Beta distribution is a conjugate distribution to the binomial distribution. In other words, 
if the prior probability is a Beta distribution and the likelihood is a binomial distribution, 
then the posterior probability will also be a Beta distribution. In this case, it is instructive 
to view the likelihood as an "operator" on the prior probability which produces as a result 
the posterior probability. When a binomial distribution "operates" on a beta distribution 
with prior hyperparameters a and {3, the result is the following posterior distribution: 
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probihypothesis; data, I) = Betaik + a,n — k + (3) (7) 
Thus, in our case where a = (3 =1: 

probihypothesis; data, I) = Betaik + 1, n — k + 1) (8) 
Thus, our posterior probability distributions for each sample are given by: 

prob{e UnX ; k u ni) = Beta(e binl ; k x + l,rai - h + 1) (9) 

prob(e bin2 \k 2 ,n 2 ) = Beta(e Un2 \k 2 + l,n 2 - k 2 + 1) (10) 
The posterior probability density for 5 is then given by: 



prob(5;k,n)= / Beta(e bin \ki + 1,^ - k x + l)Beta(e bin - 5\k 2 + l,n 2 - k 2 + 1) de bin (11) 

J — oo 

We used Monte Carlo methods in the R programming language to evaluate this integral. 
10 4 random deviates were taken from each posterior probability Beta distribution and the 
posterior probability density for 5 was determined from these. Two representative posterior 
probability densities (for the case where both samples likely share the same binomial distri- 
bution and also the case where binomial distributions differ between samples) are presented 
in Fig. O Here we present the l-a (68%) and 2-cr (95%) confidence intervals from the pos- 
terior probability density for 5 as a counterpart to the Fisher exact test likelihoods. These 
confidence intervals quantify the probable relationship between the true binary fractions of 
the two samples. For example, for the second comparison presented in Fig. [6J a sample with 
binaries detected out of 25 objects compared with a sample of 6 binaries detected out of 
23 objects, at the l-a level e bin2 for the 2nd distribution lays between e bin \ + 0.15 and e bin i 
+ 0.33. 



7.3. Brown Dwarf Binary Fraction as a Function of Mass 



We compare our measured binary fraction to tha t of more mas s ive br own dwarfs and very 
low mass stars in the Upper Sco embedded cluster. iKraus et al.l (120051 ) surveyed 12 brown 
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dwarfs and very low mass stars with ACS on HST. These objects have estim ated masses 
of 0.04 - 0.1 M and thus comprise a higher mass sample than our survey. iKraus et al. 



(120051 ) discovered three binaries in this sample, one of which (USco-109 AB) is below the 
sensitivity of our survey to detect, with a projected separation of only ~5 AU. Thus, for the 
purpo ses of comparison, we adopt a bi nary fraction of 2/1 2 = 17+g 5 % for the IKraus et al. 



(120051 ) sample. The likelihood that the IKraus et all ( 120051 ) sample is drawn from the same 
distribution as ours is 0.15, with a 1— ex B ayesian confidence interval of e bin i — e Un2 = -0.28, 
0.15. Thus, as noted by previous authors IKraus et all ( 120051 ). the binary fraction in Upper 
Sco continues to decrease with decreasing primary mass. 

This comparison is limited by the relatively small number of objects observed in Upper 
Sco. Thus, to improve statistics, we have compiled a larger list using objects from similar 
surveys of other young, nearby regions - specifically Taurus (<1 Myr, 145 pc, objects from 



Kraus et al 



20061: iKonopacky et all 120071 ) and Chamaeleon (<3 Myr, 160 pc, objects from 



Ahmic et al.l 120071 ). We include only companions that would have been detected at the 
sensitivity level of our survey and initially limit this analysis to nearby clusters (<200 pc) 
since more distant clusters (e.g. NGC 1333, IC 348, Serpens) are more than 250 pc distant 
and do not reach comparable sensitivity levels at 10 AU. All selected surveys have similar 
sensitivity levels (complete to q~0.8 at 10 AU, complete to q~0.2 - 0.3 at >20 AU) so it is 
unlikely that our survey would have discovered a binary at a separation >10 AU missed by 
these other surveys, and vice versa. We adopt 3 mass bins for this analysis: (1) high mass 
(0.07 - 0.1 M ), with 6 binaries detected out of 23 objects surveyed (6 objects from Ahmic 
et al. 2007, 5 from Kraus et al. 2006, 6 from Kraus et al. 2005, 4 from Konopacky et al. 
2007, and the two objects from the Slesnick et al. 2008 sample from the current survey), 
(2) medium mass (0.04 - 0.07 M ), with binary detected out of 18 objects surveyed 
(4 objects from Ahmic et al. 2007, 6 from Kraus et al. 2005, and 8 from Kraus et al. 
2006), and (3) low mass (<0.04 M ), with binaries detected out of 25 objects surveyed 
(7 objects from Kraus et al. 2006 and the 18 objects from the Lodieu et al. 2008 sample 
surveyed herein). We note that while a number of ad ditional bi naries are known in this 



mass range, e.g. 2 MASS 1207AB (IChauvin et al. 



2006 



Javawardhana fc Ivanovl 120061 ; 



Allers et al. 



2007 



2005). 2M 1622 (Allers 120061 : lAllers et al. 



Close et al.ll2007f ). UScoCTIO 108 



feeiar et al.ll2008h . and 2 MASS 0441 ( ITodorov et al.lboiol T survey statistics are not available 



for these objects and thus we cannot include them in our sample. Binary fractions and 
likelihoods between bins as a function of mass are presented in Table [6j As expected, the 
binary fraction decreases monotonically with primary mass. The likelihood that the low 
mass bin (<0.04 M ) objects share the same binary fraction as the high mass bin (>0.07 
M ) is less than 0.02, with a 1 — a Bayesian confidence interval of £&i n i — £un2 = 

-0.34, -0.15. 
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7.4. Brown Dwarf Binary Fraction as a Function of Age 



By ages of 1 Gyr, most of our survey objects will have cooled to become T dwarfs. 
Thus, it is interesting to compare the primordial binary fraction of these objects to the binary 
fraction of similar objects in the field. Our survey is only sensitive to companions at projected 
separations of >10 AU, however, this is a particularly interesti ng separation space to probe , 
as older field T dwarf binaries rarely have separations this large (IBurgasser et al.ll2003l . 120061 ) . 



In fact, of the 32 T dwarfs surveyed in IBurgasser et al.l (12003 



20061 ). no companions were 



detected with separation >10 AU (down to q>0.4, i.e. comparable sensitivity limits to our 
survey). This places an upper limit on the binary fraction >10 AU of 5%. Again using 
the Fischer exact test method, we found a likelihood of 1 with a very tight 1—a Bayesian 
confidence interval of euni — £un2 = -0.02, 0.07 - i.e., given the small sizes of both of these 
samples, they are very likely drawn from the same parent sample. Thus, the very low mass 
brown dwarf binary fraction appears to be similar for both young and field objects. Binary 
fractions, likelihoods, and Bayesian confidence intervals between bins as a function of age 
are presented in Table 

Do higher mass objects (>0.07 Ms U n) in young clusters also have a similar binary fraction 
(>10 AU) as their counterparts in the field? We compare the binary fraction (>10 AU 
separation) for 6 binaries discovered out of 23 young objects ( the "high mass" bin from 



the previous section) drawn from binarity surveys of Upp er Sco (IKraus et a 



( IKraus et al.ll2006t iKonopacky et al.l 120071 ). Chamaeleon f lAhmic et al.l 120071 ) and this work 



20051 ). Taurus 



with that of 1 binary (>10 AU separation) discovered from 39 field M8-L0.5 objects from 
Close et al. 2003. These two samples share a similar mass range (primary mass between 
0.07 - 0.1 M Q ), but very different wide binary fractions: 26^% for the young sample vs. 
2.61qq 6 % for the old field sample. Using the Fischer exact test, the likelihood that these two 
samples are drawn from the same binomial distribution is 0.01 with 1—a Bayesian confidence 
interval of eun\ — £un2 = 0.14, 0.32. Thus, for objects with mass > 0.07 M , there is an 
overabundance of 10-50 AU separation very low mass binaries in young clusters relative to 
the field. 

Upper Sco is a somewhat older and higher density OB association, while Taurus and 
Chamaeleon are more diffuse, younger T clusters. Thus, we also compare binary fraction 
between these two different ages and environments. Combining the sample described in the 
previous section and separating by region, we find 3 binaries detected from 34 objects in 
Taurus and Chamaeleon and 3 binaries detected from 32 objects in Upper Sco. The binary 
fraction is nearly the same between the two, although it is important to note that the sample 
in Taurus and Chamaeleon has systematically higher primary masses than that in Upper Sco 
(dominated by the 18 very low mass brown dwarfs surveyed in this paper.) Thus, given the 
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trend in binary fraction with mass, the binary fraction in Upper Sco may be considerably 
higher than in Taurus. 



7.5. Trends in Very Wide Binarity (30 — 500 AU) as a Function of Age and 

Mass 



We initially limited our statistical analysis to nearby clusters (<200 pc) since more 
distant clusters (e.g. NGC 1333, IC 348, Serpens) are more than 250 pc distant and do not 
reach comparable sensitivity levels at separations of 10-30 AU. However, including results 
from surveys of these more dis tant clusters significan tly boosts sample size. In particular, 
including the results from the iLuhman et al.l ( 120051 ) HST survey of IC 348 (2 Myr, 315 
pc) introduces 31 additional <0.1 M© objects into this analy sis. In order to match the 



achieved contrast and physical resolution of the ILuhman et al.l (120051 ) survey with those of 
nearer regions (Upper Sco, Taurus, Chamaeleon) we only consider results for separations 
from 30-500 AU. 



Including objects from lLuhman et al.l (120051 ) in our 3 mass bins from earlier sections, we 
now find: (1) in the high mass bin (0.07 - 0.1 M Q ), 3 binaries are detected out of 43 objects 
surveyed (USco-55 and USco-66 from Kraus et al. 2005 and UScol609 from this work have 
separations <30 AU and thus would not be detected at the combined sensitivity limits for 
our composite survey), (2) in the medium mass bin (0.04 - 0.07 M ), binaries detected 
out of 20 objects surveyed, and (3) in the low mass bin (<0.04 M ), binaries detected out 
of 36 objects surveyed. 

As before, the lowest mass cluster bin possesses a very similar upper limit on binarity 
as the field T dwarf bin and wide binaries seem to be rare in both the medium and low mass 
cluster bins. Comparing the high mass bin with the low mass bin, we find a likelihood of 0.25 
that these two samples are drawn from the same binomial distribution, with a 1— o Bayesian 
confidence interval of e^ n i — tuni = -0.11, -0.02. We also compare the high mass bin with 
the Close et al. 2003 sample (adjusting contrast levels, we now find binaries imaged with 
separations >30 AU out of 39 surveyed objects). According to the Fisher Exact Test, the 
likelihood that these two samples were drawn from the same binomial distribution is 0.24 (as 
opposed to 0.01 for the same bin in the nearby sample.) Thus, a significant overdensity of 
young binaries relative to the field is apparent in this sample only at moderate separations 
(10-30 AU), and not at wide (30-500 AU) separations. 
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7.6. Stability of 10-50 AU Separation Binaries in Young Nearby Starforming 

Regions 

Up to ~25% of very low mass (henceforth VLM) stars and substellar objects in young 
star forming regions may have companions at separations >10 AU. However, very low mass 
star / brown dwarf binaries with separations >15 AU are rare in the field. Of ~100 VLM 
binaries compiled at vlmbinaries.org, only ~10% have separations >15 AU. Assuming a 
binary fraction of ^ 10%, this means t hat less than 1% of field VLMs have companions at 
separations >15 AU Idose et all ( 120071 ). 



Close et all (120071 ) suggest that very wide (>50 AU) young brown dwarf binaries are 



disrupted within the first 10 Myr of their existence by interactions with stars in their natal 
cluster. To set limits on the survival time of a young wide binary in its natal cluste r , they 
adopt the analytic solution of Fokker-Plank (FP) coefficients from IWeinberg et all (119871 ) 
which describes the advective diffusion of a binary due to stellar encounters. Namely, from 
this solution, the time t* necessary to evaporate a binary with initial semimajor axis ao is: 



U ~ 3.6 x 10 £ 



r)( 



M, 



°*)( A 



\-2/ 



V, 



rel 



0.05pc" 3 A M Q n M Q J K 20fcms~ 1 ; v AU 



t)( 



do 



v-1 



(12) 



where n* is the number density of st e llar p erturbers of mass M* and relative velocity 
V re i- Using this relationship, I Close et all (120071 ) determine that young wide VLM binaries 
such as 2M 1207-39AB will not survive 10 Myr of interactions with 0.7 M stellar perturbers 
with a number density n* of 1000 pc~ 3 . Thus, I Close et all (120071 ) show that most of these 
binaries will not survive to join the field if born in a dense starforming region. W e determine 
here w hether the same is true for moderately wide 10 - 50 AU binaries. While IClose et al. 
( 120071 ) assume a number density of nearby stars of 1000 pc~ 3 , which is appropriate near dense 
core regions, it is probably too high for objects in diffuse areas of Taurus or Chamaeleon. 
Assuming a typical density of 100 pc~ 3 for for Taurus, Ophiuchus, and Upper Sco, we repeat 
this calculation for the six binaries that fall i nto our highest mass bin (specifically, CFHT-Tau 
7, C FHT-Tau 17 and C FHT-Tau 18 from iKonopackv et all (l2007h . USco-55 and USco-66 
from lKraus et all (120051 ) . and SCH 1609AB, the newly discovered binary presented herein). 
We find that all of these binaries are quite stable and will survive >10 Myr in either a 100 
pc -3 environment or a 1000 pc -3 environment (i.e. long enough to join the field population). 
An environment with stellar densities >10 4 pc -3 (equivalent to the Trapezium cluster in 
Orion) is necessary to disrupt these binaries on <10 Myr timescales. 

The existence of a significant population of these medium-separation binaries presents 
a conundrum, since very low mass stars and brown dwarf binaries with separations >15 AU 
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are rare in the field. However, the field brown dwarf populati on encompa s ses a mix of objects 
which formed in a variety of different star-forming regions. IClose et al.l ( 120071 ) suggest that 
brown dwarf binaries with separations >20 AU are found rarely in the field because they 
can only form in low-density star-forming regions, while the majority of field objects formed 
in denser initial regions where any such binary would be disrupted. However, other authors 
have suggested that most stars in the field likely fo rm in OB associations like Upper Sco 
( iKonopacky et al.l 120071 ; iPreibisch fc Mamajekll2008l ). so this is problematic. 



The existence of this population of moderately wide young brown dwarf binaries in 
lower density young clusters initially suggests that most (predominantly single) field brown 
dwarfs must form in high stellar density regions which disrupt such wide binaries by late ages. 
However, this supposition relies on our ability to distinguish between "typical" vs. "atypical" 
star-forming regions, as well as to truly disentagle the primordial vs. evolved populations. 
In other words, the evolved population is the outcome of the formation mechanism plus 
any subsequent evolution in the cluster. Different combinations of formation and subsequent 
evolution may form the same evolved population. Here, we have placed constraints on binary 
evolution within a relatively diffuse cluster environment; placing constraints on formation 
mechanism is more difficult. 

Forming brown dwarfs at all has always been a tricky prospect theoretically. Brown 
dwarf formation theories require either: (1) a mechanism to produce very low Jeans masses 



Padoan & Nordlund 


2004; 


Bate 


2009; 


Bonnell et al. 


2008 



Price fc Batd 120081 ) or (2) a method to circumvent the need for very low Jeans masses, 
(e.g. ejection, or gravitat i onal instability followed by bina ry disruption iReipurth fc Clarke 
200 ll ; IStassun et al.ll2007l ; IStamatellos fc Whitworthll2009l ) Unfortunately, more information 
is needed regarding the physical properties of these star-forming regions (i.e. measurement 
of turbulent motions, magnetic fields) to distinguish between these models. For instance, 
widespread filamentary structure has recently been observed by Herschel in very young star- 
forming clouds in Aquila and Polaris ( jMen'shchikov et al.ll2010l ). However, it is not currently 
clear what causes these filaments; if turbulence or magnetic fields are the dominant cause, 
this has significant ramifications for subsequent brown dwarf formation in these regions. 

Likely a mix of formation mechanisms are at play in any given region, the detailed 
physics of which may vary from region to region. Disentangling these physics is a difficult 
prospect and requires more information than just binary fraction. While we can rule out 
pure ejection (without any dissipation from e.g. a circumstellar disk) from our measured 
binary fraction and the existence of a significant population of >10 AU separation binaries, 
other models may produce a significant wide binary population which may be disrupted by 
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late ages in dense clusters. 

8. Conclusions 

We searched for binary companions to 20 brown dwarfs in Upper Scorpius (145 pc, 5 
Myr, nearest OB association) with the laser guide star adaptive optics system and the facility 
infrared camera NIRC2 on the 10 m Keck II telescope. This survey is the most extensive to 
date for companions to very young (5 Myr), very low mass (<40 Mj up ) cluster brown dwarfs. 
We discovered a close companion (0.14", 20.9±0.4 AU) to the very low mass object SCH 
J 1609 1837-20073523. From spectral deconvolution of integrated-light near-IR spectroscopy 
of SCH1609-2007 using the SpeX spectrograph (Rayner et al. 2003), we estimate primary 
and secondary spectral types of M6±0.5 and M7±1.0, corresponding to masses of 79±17 
Mj up and 55±25 Mj up at an age of 5 Myr and masses of 84±15 Mj up and 60±25 Mj up at an 
age of 10 Myr. 

For our survey objects with spectral types later than M8, we find an upper limit on 
binary fraction of <9% (1-0") at separations greater than 10 AU. As expected from similar 
mass binaries in the field, we find that the binary fraction (10 - 500 AU separations) appears 
to decrease monotonically with mass for young brown dwarfs. However, while proto-T-dwarfs 
(M<40 Mj up ) have a similar wide (10 - 500 AU) binary fraction as field T dwarfs, there exists 
an anomalous population of wide higher mass binaries (0.07 - 0.1 M Q primaries, separations 
of 10-50 AU) at young ages relative to older ages. 

The data presented herein were obtained at the W.M. Keck Observatory, which is oper- 
ated as a scientific partnership among the California Institute of Technology, the University 
of California and the National Aeronautics and Space Administration. The Observatory 
was made possible by the generous financial support of the W.M. Keck Foundation. The 
authors wish to recognize and acknowledge the very significant cultural role and reverence 
that the summit of Mauna Kea has always had within the indigenous Hawaiian community. 
We are most fortunate to have the opportunity to conduct observations from this moun- 
tain. B.B. was supported by Hubble Fellowship grant HST-HF-01204.01-A awarded by the 
Space Telescope Science Institute, which is operated by AURA for NASA, under contract 
NAS 5-26555. B.B. would like to acknowledge Geoffrey Mathews and Derek Kopon for help 
with observations and Adam Kraus and Eric Mamajek for useful discussions. We thank the 
anonymous referee for useful suggestions which helped improve this work. 
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Fig. 1— Left: J, H, and K s -hand images of SCH 16091837-20073523AB obtained with 
NIRC2 and the LGS AO system of the 10m Keck II telescope. North is up, east is left. Note 
that primary and companion both appear slightly elongated in the direction towards the tip- 
tilt star. The confirmed companion is at 0.144±0.002" separation and PA=15.87±0.13° with 
flux ratios of A J = 0.51±0.09, AH = 0.51±0.03, and AK S = 0.46±0.01 mag. We estimate 
masses of 47.4±11.7 Mj up and 33.5±6.0 Mj up for primary and companion respectively. 
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Fig. 2. — Left: the JHK$ colors of SCH 1609 AB compared to Upper Sco objects with M 
and L spectral types (Slesnick et al. 2008, Lodieu et al. 2008) and young field brown dwarfs 
from Cruz et al. 2009. SCH 1609ABs' colors are plotted as a red circles and are consistent 
with those of a mid to late M dwarf. Errors on SCH 1609AB photometry are shown in the 
top left corner. The DUSTY 5 Myr isochrone is plotted as a solid line. DUSTY models 
predict considerably bluer colors at these ages than is observed. Right: J-K$ vs. for 
SCH 1609AB and the same set of comparison objects. SCH 1609AB are plotted as red 
circles; combined photometry for the system is plotted as a red triangle. The DUSTY 5 Myr 
isochrone is again plotted as a solid line; while K$ band magnitudes agree with DUSTY 
predictions, colors are considerably redder than the predictions for M dwarfs. 
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Fig. 3. — Composite near-IR spectrum of SCH1609AB (black), compared to the best-fitting 
synthetic composite spectrum (red). The synthetic composite spectrum is the combination 
of UScoCTIO 75 (M6, Ardila et al. 2000; Preibisch et al. 2002) and DENIS-P J155605.0- 
210646 (M7; Martin et al. 2004, Slesnick et al. 2008). 
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Fig. 4. — Left: 5a contrast curves for 18 survey objects from lLodieu et al.l (120081 ). Noise 
levels after data reduction were calculated as a function of radius by calculating the standard 
deviation in an annulus (with width equal to approximately the FWHM of the PSF) centered 
on that radius. Noise curves were then converted to contrast in A magnitudes by dividing 
by the measured peak pixel value of the object. In general, we achieve contrasts of >4 mag 
at separations of >0.4" , sufficient to detect a 2MASS 1207 analogue at the distance of Upper 
Sco. Right: Minimum detectable absolute magnitude for the same 1 8 objects. Contrasts 
wer e converted int o abso lute magnitudes using photometry reported in lLodieu et al.l ( 2008 ) 
and lSlesnick et al.l (120081 ). and adopting a distance of 145 p c for Upper Sco. A fi lter transform 
was calculated from K to K s band using the spectra from lLodieu et al.l (120081 ) 
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Fig. 5. — Mass ratio (q) vs. separation using the DUS TY models. Con t rasts were converted 
to minimum detectable mass ratios using the models of IChabrier et al.l ( I2000I ) at an adopted 
age of 5 Myr. For the best 75% of our data, we are complete to q~0.8 at 10 AU and complete 
to q~0.2 at >20 AU. 
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p1 (0/1 8) - p2(0/32) p1 (0/25) - p2(6/23) 

Fig. 6. — Sample posterior probability distributions. Blue lines show the 1-cr confidence 
intervals on e^ni — ^Un2 = pi - p2. Left: two samples which are likely drawn from the same 
binomial distribution (specifically <0.04 M Q BDs in Upper Sco with binaries detected out 
of 18 objects, compared with field T dwarfs, with binaries detected out of 32 objects). 
The posterior probability distribution is strongly peaked at and shows little spread. Right: 
two samples which are likely drawn from different binomial distributions (specifically <0.04 
M Q Upper Sco, Taurus, and Chamaeleon objects, with binaries detected out of 25 objects, 
compared with 0.07-0.1 M Q Upper Sco, Taurus, and Chamaeleon objects, with 6 binaries 
detected out of 23 objects). The posterior probability distribution in this case peaks con- 
siderably away from 0, and is wider and flatter than the previous case. In particular, at the 
l-o" level eun2 = p2 for the 2nd distribution is between e^ni + 0.15 and euni + 0.33. 



Table 1. Known Young (<15 Myr) Very Low Mass Mass Binaries (Primary Mass <0.1 M ) 



ID 


RA 


Dec 


SpTl 


SpT2 Massl 


Mass2 


Proj. Sep. (") 


Proj. Sep. (AU) 


Ref,Notes 


Orion (400 pc, <1 Myr) 


2MASS J05352184-0546085 


05:35:21.84 


-05:46:08.5 


M6.5 


M6.5 55 M Jup 


35 M Jup 




0.04 AU 


a,b 


Taurus (140 pc, <1 Myr) 


V410-Xray3 


04:15:01. 9 C 


28:18:48. c 


M6 


M7.7 0.093 Mq 


0.047 M 


0.044±0.002" 


-6 AU 


d,e 


MHO-Tau-8 


04:33:01.1 


24:21:11.0 


M6 


M6.6 0.097 M 


0.073 M@ 


0.044±0.008" 


-6 AU 


d,f 


2MASS J04414489+2301513AB 


04:41:44.89 


23:01:51.3 


M8.5 


-20 M Jup 


5-10 M Jup 


1.105" 


15 AU 


g 


CFHT-Tau 18 


04:29:21.65 


27:01:25.95 


M6.0 


0.1 Mq 


0.06 Mq 


0.216±0.002" 


30.2 AU 


h.i 


CFHT-Tau 7 


04:32:17.86 


24:22:14.98 


M6.5 


0.07 Mq 


0.06 Mq 


0.224±0.002" 


31.4 AU 


h.i 


CFHT-Tau 17 


04:40:01.74 


25:56:29.23 


M5.75 


0.1 Mq 


0.06 Mq 


0.575±0.002" 


80.5 AU 


h.i 


FU Tau AB 


04:23:35.39 


25:03:03.05 


M7.25 


M9.25 -0.05 M Q 


-0.015 Mq 


5.7" 


800 AU 


j 


Ophiuchus (125 pc, <1 Myr) 


Oph 16AB 


16:23:36.09 


-24:02:20.9 


M5±3 


M5.5±3 -100 M Jup 


-73 M Jup 


1.7" 


212±43 AU 


k,l 


Oph 11AB 


16:22:25.21 


-24:05:13.94 


M9±0.5 


M9.5±0.5 17±* M Jup 


14±« 


1.9" 


243±55 AU 


k,l,m,n 


LkHa233 Group (325^5Qpc, ~1 Myr) 


2MASS J22344161+4041387 


22:34:41.61 


40:41:38.7 


M6 


M6 -0.1 M 


-0.1 Mq 


0.16" 


51 AU 


CO ° 








Chamaeleon (160 pc, <3 Myr) 










Cha Ha 8 


11:07:47.8 


-77:40:08 


M6.5 


0.07- 0.1 Mq 


30-35 M Jup 




1 AU 


p,q,r 


2MASS J11011926-7732383AB 


11:01:19.22 


77:32:38.60 


M7.25±0.25 M8.25±0.25 0.05 M 


0.025 M 


1.44" 


240 AU 


s 


Upper Sco (145 pc, —5 Myr) 


USco-109AB 


16:01:19.10 


-23:06:38.6 


M6 


M7.5 0.07 Mq 


0.04 Mq 


0.034±0.02" 


-5 AU 


t,u 


USco-66AB 


16:01:49.66 


-23:51:07.4 


M6.0 


M6.0 0.07±0.02 M Q 


0.07±0.02 Mq 


0.07" 


10.19±0.07 AU 


t.u 


USco-55AB 


16:02:45.60 


-23:04:49.8 


M5.5 


M6.0 0.10±0.03 M Q 


0.07±0.02 M 


0.12" 


17.63±0.09 AU 


t.u 


UScoCTIO108AB 


16:05:53.94 


-18:18:42.7 


M7 


M9.5 60±20 M Jup 


14±| M Jup 


4.6±0.1" 


-670 


v,u 


R Corona Australis (-130 pc, -0.5-10 Myr) 


DENIS-P J185950.9-370632 


18:59:50.9 


-37:06:32 


M8±0.5 


0.017 Mq 


0.013 Mq 


0.06" 


7.8 AU 


w 


TW Hydra (-30-70 pc, -12 Myr) 


2MASS J1207334393254 


12:07:33.40 


39:32:54.0 


M8 


L5-L9.5 -25 M Jup 


5±2 M Jup 


0.78" 


55 AU 





HStassun etaD l|2006h 
HStassun et al. 
c epoch 1950 coordinates 
^Kraus et all bOOrj'l 
fetrom fc Stroml l|l994h 



iBriceno et~aH l|l998h 
^Todorov et alj faOld) 
' iKonopackv et al, | j2007h 
iGuieu et all <2006h 
jLuhman et al.l j2009h 
HAllerj <2006h 
IClose et al] ll2007h 
" faavawardhana fc Ivanovl (2006) 
"estimated age~5 Myr 
* lAlIers et all teOOSh 
r jjoergensl J2006h 
' jjoergens fc Mulled l l2007f) 
" Ijoergens et al. I <201Ch 
' iLuhmanl J2004) 
iKraus et al] ||2005|) 
HArdila et all hOOdt) 
HBeiar et all ||2008|) 
HBouv et all j2004) 
HChauvin etal] l|2005h 
HMohantv et al.l fcOOTh 



Co 
to 
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Table 2. 


Objects Observed 










ID 


Right Ascension 


Declination 


SpT 


.1 


H 


K 


flcCOsS^ 




Objects from Lodieu et al. 2008 sample 


USco J155419.99-213543.1 


15:54:19.99 


-21:35:43.1 


M8 


14.93 


14.28 


13.71 


-14 


-18 


USco J160603.75-221930.0 


16:06:03.75 


-22:19:30.0 


L2 


15.85 


15.10 


14.44 






USco J160606.29-233513.3 


16:06:06.29 


-23:35:13.3 


L0 


16.20 


15.54 


14.97 




- 4 


USco J160714.79-232101.2 


16:07:14.79 


-23:21:01.2 


L0 


16.56 


15.83 


15.07 




- 4 


USco J160723.82-221102.0 


16:07:23.82 


-22:11:02.0 


LI 


15.20 


14.56 


14.01 


-11 


-31 


USco J160727.82-223904.0 


16:07:27.82 


-22:39:04.0 


LI 


16.81 


16.09 


15.39 






USco J160737.99-224247.0 


16:07:37.99 


-22:42:47.0 


L0 


16.76 


16.00 


15.33 






USco J160818.43-223225.0 


16:08:18.43 


-22:32:25.0 


L0 


16.01 


15.44 


14.70 






USco J160828.47-231510.4 


16:08:28.47 


-23:15:10.4 


LI 


15.45 


14.78 


14.16 


-12 


-13 


USco J160830.49-233511.0 


16:08:30.49 


-23:35:11.0 


M9 


14.88 


14.29 


13.76 


-5 


-12 


USco J160843.44-224516.0 


16:08:43.44 


-22:45:16.0 


LI 


18.58 


17.22 


16.26 




- 12 


USco J160847.44-223547.9 


16:08:47.44 


-22:35:47.9 


M9 


15.69 


15.09 


14.53 





-20 


USco J160918.69-222923.7 


16:09:18.69 


-22:29:23.7 


LI 


18.08 


17.06 


16.16 




- 8 


USco J161047.13-223949.4 


16:10:47.13 


-22:39:49.4 


M9 


15.26 


14.57 


14.01 


-15 


-24 


USco J161228.95-215936.1 


16:12:28.95 


-21:59:36.1 


LI 


16.41 


15.56 


14.79 






USco J161302.32-212428.4 


16:13:02.32 


-21:24:28.4 


L0 


17.17 


16.37 


15.65 






USco J161441. 68-235105.9 


16:14:41.68 


-23:51:05.9 


LI 


16.07 


15.34 


14.62 






USco J163919.15-253409.9 


16:39:19.15 


-25:34:09.9 


LI 


17.20 


16.39 


15.61 


-1 


-17 


Additional Objects 


SCH J16091837-20073523 


16:09:18.37 


-20:07:35.23 


M7.5 


13.00 


12.37 


12.01 






SCH J16224384-19510575 


16:22:43.84 


-19:51:05.75 


M8 


12.35 


11.61 


11.15 







HLodieu et ail \200l\) 
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Table 3. Observations 



ID 


Observation Date 


Filter 


Exposure Time 


Median Strehl 


Median FWHM 


Objects from Lodieu et al. 2008 sample 


USco J155419. 99-213543.1 


2009-05-30 


K s 


7x60 s 


0.10 


99 mas 


USco J160603. 75-221930.0 


2007-07-17 


Ks 


11x15 s 


0.19 


66 mas 






J 


9x30 s 


0.03 


66 mas 






11 


9x15 s 


0.06 


57 mas 




2008-07-27 


K s 


11x15 s 


0.31 


55 mas 


USco J160606. 29-233513. 3 


2009-05-29 


K S 


6x60 s 


0.15 


81 mas 


USco J160714.79-232101.2 


2009-05-29 


K S 


6x60 s 


0.10 


80 mas 






J 


6x60 s 


0.02 


95 mas 






H 


6x60 s 


0.05 


92 mas 


USco J160723. 82-221102.0 


2008-07-27 


K S 


14x15 s 


0.25 


67 mas 






J 


12x30 s 


0.03 


82 mas 






11 


12x15 s 


0.11 


66 mas 




2009-05-30 


K S 


7x60 s 


0.13 


82 mas 


USco J160727.82-223904.0 


2008-07-27 


K S 


11x15 s 


0.21 


62 mas 






J 


9x30 s 


0.02 


74 mas 


USco J160737.99-224247.0 


2009-05-29 


K S 


6x60 s 


0.11 


96 mas 


USco J160818. 43-223225.0 


2009-05-29 


K S 


6x60 s 


0.24 


66 mas 


USco J160828. 47-231510. 4 


2007-07-17 


K S 


11x15 s 


0.13 


82 mas 


USco J160830. 49-233511.0 


2009-05-30 


K S 


7x60 s 


0.06 


130 mas 


USco J160843. 44-224516.0 


2009-05-29 


K S 


6x60 s 


0.16 


78 mas 






J 


6x60 s 


0.02 


70 mas 


USco J160847.44-223547.9 


2007-07-17 


K S 


11x15 s 


0.23 


68 mas 


USco J160918. 69-222923. 7 


2008-07-27 


K S 


11x15 s 


0.13 


67 mas 


USco J161047.13-223949.4 


2007-07-17 


K S 


11x15 s 


0.19 


73 mas 


USco J161228. 95-215936.1 


2008-07-27 


K S 


11x15 s 


0.16 


67 mas 


USco J161302. 32-212428. 4 


2009-05-30 


K S 


7x60 s 


0.16 


79 mas 


USco J161441. 68-235105. 9 


2008-07-27 


K S 


11x15 s 


0.16 


67 mas 


USco J163919. 15-253409. 9 


2008-07-27 


K S 


11x15 s 


0.19 


62 mas 






J 


9x30 s 


0.02 


69 mas 


Additional Objects 


SCH J16091837-20073523 


2009-06-30 


K S 


6x20 s 


0.11 


75 mas 






J 


6x20 s 


0.01 


80 mas 






11 


6x20 s 


0.04 


71 mas 


SCH J16224384-19510575 


2009-05-30 


K S 


6x10 s 


0.29 


59 mas 
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Table 4. Properties of the SCH1609-2007AB System 



Primary Secondary 



Distance 
Age 

Separation 
Position Angle 
A J (mag) 
AH (mag) 
AK S (mag) 
J (mag) 
H (mag) 
K s (mag) 
J - Ks (mag) 
J — H (mag) 
H - K s (mag) 

Spectral Type 

T eff 

Estimated Mass (5 Myr) 
Estimated Mass (10 Myr) 



145±2 pc a 
5 Myr b 

0.144±0.002" (20.9±0.4 AU) 
15.87±0.13° 

0.51±0.09 
0.51±0.03 
0.46±0.01 
13.53±0.09 c 14.04±0.09 
12.90±0.04 c 13.41±0.04 
12.56±0.03 c 13.01±0.03 
0.97±0.09 1.03±0.09 
0.63±0.10 0.63±0.10 
0.34±0.05 0.40±0.05 
-2.04±0.12 -2.23±0.12 

M7±0.5 M6±1.0 
2990±60 K 2850±170 K 
79±17 M Jup 55±25 M Jup 
84±15 M Jup 60±25 M Jup 



Preibisch et all (|2002l ) 



IPreibisch fc Zinneckerl (|l999h 
c from 2MASS 
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Table 5. Measured K s Contrast and Minimum Detectable Mass Ratios 



ID 


Observation Date 


Amag(0.07 ) 


q(0.07 ) 


Amag(0.2 ) 


q(0.2") 


Amag(0.5") 


q(0.5") 


USro T1 55419 99-213543 1 


2009-05-30 


0.41 


0.87 


3.67 


0.24 


5.54 


0.11 


USro T1 60603 75-221930 


2008-07-27 


1.83 


0.46 


5.01 


0.13 


5.88 


0.08 


USro T1 60603 75-221930 (1 


2007-07-1 7 


1.17 


0.59 


4.40 


0.16 


5.42 


0.10 


USco J160606. 29-233513. 3 


2009-05-29 


0.65 


0.71 


3.79 


0.20 


4.84 


0.13 


USco J160714. 79-232101. 2 


2009-05-29 


0.70 


0.71 


3.89 


0.19 


4.66 


0.14 


USco J160723.82-221102.0 


2008-07-27 


1.23 


0.64 


4.95 


0.14 


5.98 


0.08 


USco J160723.82-221102.0 


2009-05-30 


0.62 


0.80 


3.74 


0.23 


5.53 


0.11 


USco J160727.82-223904.0 


2008-07-27 


1.36 


0.54 


4.28 


0.16 


4.55 


0.14 


USco J160737.99-224247.0 


2009-05-29 


0.42 


0.86 


3.45 


0.23 


4.22 


0.16 


USco J160818. 43-223225.0 


2009-05-29 


1.16 


0.58 


4.53 


0.15 


5.61 


0.08 


USco J160828.47-231510.4 


2007-07-17 


0.59 


0.80 


4.14 


0.19 


5.26 


0.12 


USco J160830.49-233511.0 


2009-05-30 


0.14 


0.95 


2.91 


0.32 


5.11 


0.13 


USco J160843.44-224516.0 


2009-05-29 


0.78 


0.72 


3.20 


0.27 


3.47 


0.23 


USco J160847.44-223547.9 


2007-07-17 


1.23 


0.56 


4.68 


0.14 


5.51 


0.10 


USco J160918.69-222923.7 


2008-07-27 


1.11 


0.63 


3.46 


0.23 


3.66 


0.21 


USco J161047. 13-223949.4 


2007-07-17 


0.82 


0.75 


4.49 


0.17 


5.66 


0.10 


USco J161228. 95-215936.1 


2008-07-27 


1.08 


0.60 


4.23 


0.17 


4.78 


0.13 


USco J161302. 32-212428.4 


2009-05-30 


0.69 


0.76 


3.91 


0.18 


4.70 


0.11 


USco J161441. 68-235105.9 


2008-07-27 


1.48 


0.52 


4.52 


0.15 


5.18 


0.11 


USco J163919. 15-253409.9 


2008-07-27 


1.38 


0.56 


3.98 


0.18 


4.36 


0.14 



Table 6. Statistical Sample Comparison as a Function of Mass 



Sample 1 


Sample 2 


Likelihood l-cr CI on 5 = £t,; n i — £(,i n 2 


2-cr CI on 5 = e binl — e bin2 




Upper S< 


:o Mass Comparison (10 


- 1000 AU scparat 


ions) 


<0.04 M 

/ 18, <9% 


0.04-0.1 M 
2 / 12, 17±g 5 % 


0.15 


-0.28, -0.15 


-0.42, 0.04 


Upper Sco, Taurus, and Chamaeleon Mass Comparison (10 - 1000 AU separations) 


<0.04 M 

/ 25, <7% 


0.07-0.1 M 
6 / 23, 26±"% 


0.01 


-0.34, -0.15 


-0.44, -0.07 


<0.04 M 
/ 25, <7% 


0.04-0.07 M 
/ 18, <9% 


1.0 


-0.06, 0.04 


-0.15, 0.10 


0.04-0.07 Mq 

/ 18, <9% 


0.07-0.1 M 
6 / 23, 26l^% 


0.03 


-0.33, -0.13 


-0.43, -0.04 


Upper S 


co, Taurus, Chamaeleon, and IC 348 Mass Comparison (30 - 


1000 AU separations) 


<0.04 M 
/ 36, <4.8% 


0.07-0.1 M 
3 / 43, 7±l% 


0.25 


-0.11, -0.02 


-0.17, 0.03 


<0.04 M 
/ 36, <4.8% 


0.04-0.07 M 

/ 20, <8% 


1.0 


-0.06, 0.02 


-0.14, 0.07 


0.04-0.07 M 
/ 20, <8% 


0.07-0.1 M 
3 / 43, 7±l% 


0.54 


-0.10, 0.01 


-0.16, 0.09 
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Table 7. Statistical Sample Comparison as a Function of Age 



Sample 1 Sample 2 Likelihood 1-a CI on <5 = tn n i — tbinl 2-<r 


CI on 5 — e binl — e bin 2 


Upper Sco vs. Field (10 - 1000 AU separations) 


<0.04 Mq BDs in Upper Sco Field T Dwarfs 

/ 18, <9% / 32, <5% 1.0 -0.02, 0.07 


-0.08, 0.15 


Upper Sco, Taurus, and Chamaeleon vs. Field (10 - 1000 AU separations) 


0.07-0.1 M Q Cluster BDs Field M and L Dwarfs 

6 / 23, 26^"% 1 / 39, 2.6 +5 i 1 % 0.01 0.14, 0.32 


0.06, 0.43 


Upper Sco, Taurus, Chamaeleon, and IC 348 vs. Field (30 - 1000 AU separations) 


0.07-0.1 M Cluster BDs Field M and L Dwarfs 

3 / 43, 7±2% / 39, <4.4% 0.24 0.02, 0.11 


-0.02, 0.17 



